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ABSTRACT 

We analyse a sample of 69 QSOs which have been randomly selected in a complete 
sample of 104 QSOs {R ^ 18, 0.142 < z < 0.198, S < 10°). 60 have been observed 
with the NTT/SUSI2 at La Silla, through two filters in the optical band (WB#655 
and V#812), and the remaining 9 are taken from archive databases. The filter V#812 
contains the redshifted H(] and forbidden [OIII] emission lines, while WB#655 covers 
a spectral region devoid of emission lines, thus measuring the QSO and stellar con- 
tinua. The contributions of the QSO and the host are separated thanks to the MCS 
deconvolution algorithm, allowing a morphological classification of the host, and the 
computation of several parameters such as the host and nucleus absolute V-magnitude, 
distance between the luminosity center of the host and the QSO, and colour of the 
host and nucleus. We define a new asymmetry coefficient, independent of any galaxy 
models and well suited for QSO host studies. The main results from this study are: 
(i) 25% of the total number of QSO hosts are spirals, 51% are ellipticals and 60% 
show signs of interaction; (ii) Highly asymmetric systems tend to have a higher gas 
ionization level (iii) Elliptical hosts contain a substantial amount of ionized gas, and 
some show off-nuclear activity. These results agree with hierarchical models merger 
driven evolution. 

Key v^rords: quasars: general, galaxies: interactions, active, fundamental parameters. 



1 INTRODUCTION 

Over the past decade, many observations and studies have 
been put forwards in order to reach a better understanding 
of the interrelations between the QSOs and their hosts. This 
has been made possible thanks to the availability of high 
resolution space based data, as well as ground based imag- 
ing. Several strong correlations have been found, such as the 
famous black hole mass - bulge stellar velocity dispersion 
relation for q uiescent galaxies ("Ferrarese & Merrit^ |2000| : 
iGebhardt et a l. 2000; Bcrnardi ct al. 2007), extended after- 
wards to active galaxies (jMcLure and Dunlod l200ll . l2002l : 
[m arconi and Huntll2003h . The most popular global picture 
that has emerged is to place the QSO phenomenon in an 
evolutionary context, in which galaxies evolve hierarchically 
by successive gravitational interactions or mergers, allow- 
ing more gas to reach the central regions and possibly trig- 
ger the QSO phase. Even if the observational basis of this 
idea remains unsecure, convincing hints of merger induced 
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QSO acti vity were given using eith er hydrodynamical sim- 
ulati ons dHopkins et all l2005l l200d), hi gh re s olution imag- 
ing llBahcall et al.lll997l : iBen nert et al. 2008': 'L etawe et all 
l2008bl) , or 2D-spect r oscop y (Letawe ct al. 200'^). On the 
other hand, ISchmitj (|200l[ ) argued, by comparing samples 
of active and non active galaxies, that the percentage of 
merging or gravitationally interacting systems is not higher 
in QSOs than in other types of active or non-active galax- 
ies. However, this study, based on low resolution images 
(1.7arcsec/pxl), does not take into account faint tails or com- 
pact fea t ures i n the host typical of recent merging activity. 
iLi et al. I ((2008), with a sample of 10^ low-redshift galaxies, 
find a strong correlation between the presence of a close 
companion and the star formation, but not between close 
companions and AGN activity. The puzzle is thus far from 
being solved. 

In this framework, it is important (1) to find peculiar cases in 
specific stages of evolution, such as the famous HE0450-2958 
ijMagain et al.l[2005i ) , which would allow to better assess the 
nature an d creation mechan i sms of the black ho le-bulge co- 
evolution (|Elbaz et al.ll2009l : Ijahnke et alll2009l) and (2) to 
extract, in well resolved samples, correlations between ob- 
servables as well as the proportion of interacting systems. 
In order to bring new insights on those issues, we study a 
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sample of 69 QSOs extracted from a complete sample of 104 
QSOs drawn from different catalogues. Section 2 explains 
the construction of the sample and the observations. Sec- 
tion 3 contains a description of how the contributions of 
the host and QSO are separated. The parameters derived 
thanks to this separation (QSO and host magnitudes, cen- 
ter of luminosity, asymmetry coefficient), the correlations 
found between them, and subsequent discussion are given 
in Sections 4 and 5. Some conclusions are drawn in Section 
6. Finally, a few particularly interesting cases are analysed 
with more scrutiny in the Appendix. 



2 SAMPLE AND OBSERVATIONS 
2.1 Sample 

We selected from t he m ain catalogue^ av ailable 

dVeron-Cettv fc VeronI [20061: [Schnc ider et~al] l2005l: 
Hewitt fc Burbidgdll99"3l: IWisotzki et al. 1 1200(1 " iGreen et al.l 



1986l l all the brightest QSOs {R < 18) with S ^ 10° and 
0.142 ^ z ^5 0.198. This redshift range was chosen to enable 
the observation of low redshift quasars through two specific 
filters, one including emission lines, the other only the QSO 
and stellar continua. This provided a sample of 104 QSOs, 
9 of which had already been observed at high resolution 
(see Table [1] for references), and 60 of which were observed 
in the context of the present study. Those 69 QSOs form 
the sample analysed here. Even if only 66% complete, we 
expect our sample to be statistically relevant and large 
enough to infer robust proportion measurments. Indeed, the 
35 remaining QSOs could not be observed only because of 
bad weather conditions during parts of the two runs, which 
is a selection independent of the QSOs characteristics, and 
thus devoid of any bias. 

2.2 Observations 

The observations of the sample was made with NTT/SUSI2, 
the Superb-Seeing Imager, a direct imaging camera opti- 
mised for periods of good seeing at the ESO/La Silla obser- 
vatory, during the run A 077.B-0229 in August 2006 and run 
B 078.B-0081 in February 2007. We observed through two 
filters, V#812 and WB#665. The complete list of targets 
is given in Table [21 The observation strategy was motivated 
by two major points in the understanding of the QSO-host 
interactions. First of all, good resolution and sampling are 
necessary for detecting the host and revealing its morphol- 
ogy. SUSI2, in its 2 * 2 binning mode, offers a sampling of 
0.161 arcsec/pxl, which translates in our redshift range into 
~ 0.45 kpc/pxl. Exposure times have been estimated with 
the NTT ETC in order to reach a S/N high enough to detect 
galaxies departing by as much as 3(j from the mean magni- 
tude relation between QSO and host. Each QSO observation 
was divided into 3 or 4 exposures in each filter in order to 
avoid saturation and to efficiently remove bad CCD pixels 
and cosmic ray hits. This set up, with a typical seeing of 
0.6, gives us observations deep enough to infer the major 
galactic morphology. Non photometric observing conditions 
were reported for 8 objects. 

Secondly, our previous study l|Letawe et al.l l2008bl ) shows 
that the distribution of ionized gas in the host does not 
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Figure 1. VLT FORSl spectrum of the host galaxy of the QSO 
HE1434 — 160 after the separation from the QSO, at a redshift 
of ^ = 0.144 llLetawe et al.ll2008ah . together with the response 
curves of the two filters used for probing the gaseous and stellar 
content. 



necesserally match the structure of the stars distribu- 
tion . For instance, the host of HE0354-5500, analysed in 
iLetawe et al] (|2008bl) contains gas ionized by the QSO even 
in remote regions devoid of stars. This motivates to study 
the gas and stellar content separately for each host. That 
is why all the 60 SUSI2 QSOs have been observed through 
two filters: V#812 which, in our redshift range, contains the 
two forbidden [OIII] A4959 and [OIII] A5007 Unes along with 
H/3 A4861. Consequently, this filter allows to map the mor- 
phology of the ionized gas in the host, as we demonstrate in 
Section [5.21 The filter WB#665 contains a region devoid of 
emission lines, in which the stellar continuum should thus be 
the only contributor, allowing to map the morphology of the 
stellar content of the host. The presence of the stellar contin- 
uum in the V#812 and its influence on the observed intensi- 
ties are discussed in Section [5.21 The filters response curves 
are plotted along with a typical QSO host spectrum in Fig.[T] 
All observations were fiatfielded and cleaned from bad pixels 
using PYRAF tools. Throughout the paper, we adopt the 
following cosmological parameters: Hq = 71km s~^Mpc~^, 
Qm = 0.27 and fix = 0.73. 



3 IMAGE ANALYSIS 

3.1 NTT/SUSI2 data processing 

3.1.1 Deconvolution method 

The whole analysis of host galaxy and QSO properties re- 
lies on a good separation of those two components. We 
use a deconvolution method based on the MCS algorithm 
(Ma gain et al.l[]~998l ). which is known to be particularly well 
suited to separate point so u rces from a diffuse background 
(|Letawe et al.ll2007l ! [2008al lbl: lMagain et al.ll2005l ). Its princi- 
ple is to produce images with an improved resolution with- 
out trying to reach an infinite one. Other techniques, which 
deconvolve an image with its true point spread function 
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Name 


RA 


Dec 


z 


Filter 


Ref. 


MviQSO) 


My (host) 


Morph. 


Interac? 


PHL909 


00 54 32 


14 29 58 


0.171 


F606W 


(1) 


-23.64 


-21.74 


Ell. 


N 


0205+024 


02 05 14.53 


02 28 42.7 


0.155 


F606W 


(1) 


-23.74 


-19.84 


Spir. 


N 


PKS 0736+017 


07 36 42.49 


01 44 00.1 


0.191 


F675W 


(2) 


-23.97 


-22.78 


Spir 


Y 


PKS 1020-103 


10 20 04.2 


-10 22 33.6 


0.197 


F675W 


(3) 


7 


-21.29 


Ell. 


N 


MC 1635+119 


16 35 25.88 


11 55 46.4 


0.146 


F606W 


(4) 


-22.74 


7 


7 


Y 


PG 2349-014 


13 49 22.3 


-01 25 54 


0.174 


F675W 


(3) 


7 


-22.93 


7 


Y 


3C273 


12 29 06.7 


02 03 08 


0.158 


F606W 


(1) 


-26.34 


-22.84 


Ell. 


N 


HE1405-1545 


14 08 24.5 


-15 59 28 


0.194 


B 


(5) 


-24.2 


-23.3 


Spir. 


Y 


HE1434-1600 


14 36 49.6 


-16 13 41 


0.144 


F606W 


(6) 


-23.82 


-22.87 


Ell. 


Y 



Ta ble 1. List o f the QSOs already observ e d and main host properties. Refe rence: (l):[B ahcaII et al.l lll"997ft , (2): lDunIop et aL I lEool), 
(31: iKim et~al] (|2008h . (4): ICanalizo et al.l |(20o3); f5l: Ijahnke et al.l l|2004ah : (6l: lLetawe et alj (|2008lJ ). Objects from (1) to (4) were 
observed with HST/ACS/WFC, (5) with NTT/EFOSC2 and (6) with HST/ACS/HRC. A question mark in MviQSO) is when both 
absolute and apparent V-magnitudes are unknown. 



(PSF) and attempt to reach arbitrarily high resolutions, 
tend to result in artefacts, deblending problems, flux errors, 
and astrometry errors in the deconvolved image. In the MCS 
technique, the user chooses some finite value for the resolu- 
tion of the deconvolved output image by defining an output 
PSF smaller than the original image's PSF. The user then 
chooses an appropriate pixel scale for the deconvolved out- 
put image in order to sample the chosen PSF at least at 
the Nyquist level. By fully sampling the chosen PSF, the 
MCS algorithm eliminates artefacts in the deconvolved im- 
ages, improves deblending of point sources, and reduces flux 
and position errors. For our purpose, the deconvolved im- 
age has a pixel size times smaller than in the original image 
and a PSF of Gaussian shape, with a Full Width at Half 
Maximum (FWHM) of 2 CCD pixels (4 in the deconvolved 
images). Another advantage of the method is that the point 
sources have a known shape in the deconvolved image and 
can be explicitely separated from the diffuse components, 
which allows to get an image of the host galaxy uncontam- 
inated by the QSO. A particularity of this method in the 
present context of QSO host analysis is that it introduces 
no model or prior assumption for the shape of the host, as 
the galaxy is represented by a numerical diffuse component. 

3.1.2 PSF construction 

A crucial point in the deconvolution process is the PSF con- 
struction: the more accurate the PSF at the position of the 
QSO on the CCD, the better the deconvolution. We con- 
struct the PSF as follows: 

• Some stars (between 1 and 4) are selected on the same 
frame as the QSO, as close to it as possible and with sim- 
ilar brightness. Indeed, a star as luminous as the QSO has 
the same S/N and possibly suffers from the same deviations 
from linearity. Moreover, a star far from the QSO is more 
likely affected by PSF variations across the fleld, especially 
if located near the border. 

• The kernel s{x) of the deconvolution is first approxi- 
mated by a Moffat profile (better reproducing the wings of 
the PSF than a Gaussian profile). Because real PSFs have 
more complex structure than this analytic function, we add 
to it a numerical component whose sole purpose is to im- 
prove the fit with real stars. This kernel (Moffat profile + 
numerical background), after convolution with the PSF of 
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Figure 2. Upper-Left: Different PSF stars are selected on the 
CCD. Upper-Right: A Moffat profile is adjusted on all the stars si- 
multaneously. Bottom-Right: Numerical background added to im- 
prove the PSF. Bottom-Left: Residuals showing only faint struc- 
tures particular to each star. 



the deconvolved image (i.e. the chosen Gaussian profile) is 
then simultaneously fitted on the selected PSF stars. A de- 
tailed descr iption of the way th e kernel is constructed can 
be found in lMagain et all (|2007[) . 

An example of the whole PSF construction process is given 
on Fig. [2] Only very faint structures, different for each star, 
are left in the residuals (differences between the model and 
the data for each pixel, divided by the standard deviation in 
this pixel). It shows that the relevant available information 
(i.e. common to all selected stars) is taken into account in 
s(f). 

3.1.3 Simultaneous deconvolution 

After the construction of the PSF for each QSO exposure, 
the simultaneous deconvolution of the different images of a 
given QSO observed through a given filter can be carried 
out. The QSO image is decomposed into its point source 
component and a numerical diffuse background identical for 
all exposures, thus containing the full available information 
on the host. 

The parameters to fit are thus: the intensity and cen- 
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Name 






RA 




Dec 




z 


Run 


V (s) 


WB (s) Morphology Interaction? 




00 


57 


09, 


.9 


14 


46 


10 


0. 


,172 




1500 


1800 


Ell. 


Y 


ni 1 1 1 fi-i ni fim r 


01 


11 


10, 


.0 


-10 


16 


32 


0. 


.179 




750 


930 


Spir. 


Y 




01 


18 


45, 


.5 


13 


33 


27 


0, 


189 




2100 


2625 


Ell. 


N 




01 


55 


30, 


.0 


-08 


57 


04 


0, 


165 




640 


800 


Und. 


Y 


091 91 ^ 07*^71 Q K 


02 


12 


18, 


.3 


-07 


37 


20 


0, 


174 




1350 


1650 


Spir. 


N 


091 '^fi0-^0nzl99iS 7 


02 


13 


59, 


.8 


00 


42 


27 


0, 


182 




1350 


1800 


Und. 


Y 


09'^0n7-U0n9'^9'^ 


09 


■^0 


07, 


.0 


no 


25 


25 


u. 


1 Q8 




1 900 


1600 


Und. 


Y 


0*^991 ^-UOO'^'^I 4 


03 


22 


13, 


.9 


00 


55 


13 


0. 


.185 


\ 


640 


800 


Ell. 


Y 




1 


1 

-LU 


44, 


.5 


uu 


43 


31 


u. 


1 78 

. J- / o 




O 1 \J 


720 


Spir. 


N 


1 1 ^7nfi-uni ^Q47 Q 


11 


37 


06, 


.8 


01 


39 


48 


0. 


.193 




1536 


1720 


Spir. 


N 


1 99'^'^A_n947'=i7 9 


12 


9^; 


34, 


.8 


09 
-uz 


47 


57 


n 

u. 






1 9R0 
IZUU 


1630 


EII. 


N 


1 fil ^s2-no97'^n 


16 


15 


32, 


.3 


00 


27 


30 


0, 


,146 


B 


1260 


1310 


Ell. 


N 


90f^0'^9 0701 "^1 9 


20 


50 


32, 


.3 


-07 


01 


31 


0, 


168 




2400 


2800 


Spir. 


N 


9*^1 71 9 nO'^fiO'^ f\ 


23 


17 


11, 


.8 


-00 


36 


04 


0, 


186 




1940 


2320 


Und. 


Y 


9*^99(^0 nO'^'^'^Q 


23 


23 


00, 


.0 


-00 


53 


59 


0, 


.150 




1260 


1440 


Spir. 


N 


9'^'^1 01 0Q1 


23 


51 


56, 


.1 


-01 


09 


13 


0, 


174 




400 


500 


Ell. 


N 


C\ 0099 9044 


00 


9"^ 


08, 


.4 


9(1 
-zu 


27 


35 


A 

u. 


1 70 


\ 


240 


300 


Ell. 


N 


CT 289 


01 


00 


39, 


.5 


-25 


38 


28 


0. 


.158 




600 


800 


Ell. 


Y 




01 


34 


42, 


.7 


-41 


36 


13 


0. 


.172 




600 


750 


Und. 


Y 


A/T^ 10^09-97^7 


10 


32 


36, 


.1 


-28 


13 


26 


0. 


.148 


B 


600 


750 


Ell. 


N 


PTC^l 1 941 -"^QQ 


12 


44 


29, 


.4 


-40 


12 


46 


0, 


.191 


B 


1816 


2160 


Und. 


Y 


MS 1 s'iQi 4-n4sn 


14 


01 


36, 


.7 


04 


16 


25 


0, 


.163 


B 


1272 


1590 


Ell. 


Y 


O 1 491 001 


14 


24 


03, 


.8 


-00 


26 


58 


0, 


151 


B 


600 


750 


Spir. 


N 




17 


28 


19, 


.9 


-14 


15 


56 


0, 


.184 


A 


60 


60 


Und. 


Y 


A/T^ 9007K '^^99 


20 


11 


08, 


.8 


-36 


13 


10 


0, 


177 




420 


520 


Ell. 


Y 


O 2240-241 1 


22 


43 


40, 


.9 


-23 


55 


16 


0. 


.184 


A 


1050 


1200 


Spir. 


N 


O 99^9-94*^4 


22 


55 


25, 


.0 


-24 


18 


30 


0. 


.147 




450 


240 


Und. 


Y 


fiOZ T2S2Q27-2Q'^8 


23 


29 


27, 


.4 


-29 


38 


47 


0. 


.193 


A 


560 


720 


Spir. 


N 


TTF, 0097- "^1 1 S 


00 


29 


37, 


.3 


-31 


02 


10 


0. 


.145 


A 


900 


1125 


Ell. 


Y 


TTF, 01 0J<-^499 


01 


10 


38, 


.4 


-54 


06 


40 


0. 


.186 


A 


750 


900 


Ell. 


N 


TTF, 01 4R-'^7'^'^ 

III J \jA.'-t\j-0 1 oo 


01 


48 


21, 


.1 


-37 


40 


20 


0, 


.147 


A 


600 


750 


Spir. 


N 


■RF, 099fi '^90Q 


02 


27 


59, 


.6 


-51 


56 


32 


0, 


145 


A 


1200 


1530 


Und. 


Y 


FTF 0997 41 9*^ 


02 


29 


13, 


.3 


-41 


10 


10 


0, 


143 


A 


1200 





Ell. 


Y 


FTF 09'^0 4400 


09 


•^9 


23, 


.0 


-^o 


47 


55 


u. 


1 fi8 
. lUO 


A 


^^40 


600 


Ell. 


Y 


HF, 0441-2896 


04 


43 


20, 


.7 


-28 


20 


52 


0. 


.155 


A 


150 


180 


Ell. 


Y 


WF 1 1 ni OQi^Q 


11 


04 


16, 


.7 


1 n 

-lU 


16 


08 


u. 


1 8R 
.-LOU 


B 


1 OfiO 
-LUUU 


1260 


Ell. 


N 


TTF 1 909-0'=i01 


12 


04 


53, 


.0 


-05 


18 


13 


0. 


.169 


B 


1332 


1668 


Und. 


Y 


TTF 1911-1 QOi^ 


12 


14 


03, 


.4 


-19 


21 


43 


0. 


.148 


B 


1590 


1935 


Spir. 


Y 


HF, i2Sfi-2nni 


12 


39 


01, 


.7 


-20 


17 


30 


0. 


.196 


B 


1248 


1500 


Ell. 


Y 


TTF 1 9 "=1-04^7 


12 


58 


31, 


.0 


-04 


53 


49 


0. 


.172 


B 


1104 


1380 


Ell. 


Y 


NF 1 9"^^ 91 "^Q 


12 


59 


02, 


.4 


-21 


55 


38 


0, 


146 


B 


1200 


1500 


Ell. 


Y 


FTF 1 "^OO OR'^7 


13 


02 


46, 


.7 


-07 


13 


55 


0, 


181 


B 


1470 


1830 


Spir. 


Y 


NF 9'^4'^-'^Q'^Q 


23 


48 


12, 


.1 


-39 


23 


07 


0, 


196 


A 


1090 


1125 


Spir. 


N 


\jkjO\j 0\)0 


00 


56 


15, 


.8 


-36 


22 


17 


0, 


.162 


A 


800 


1000 


Ell. 


Y 


01 'i9-l-077 


01 


32 


31, 


.7 


07 


43 


47 


0, 


.147 


A 


1050 


1320 


Spir. 


Y 


091 "^-484 


02 


13 


52, 


.6 


-48 


26 


55 


0. 


.168 


A 


1080 


1260 


Ell. 


Y 


0*^^74-1 07 


03 


57 


27, 


.1 


10 


46 


48 


0. 


.182 


A 


544 


680 


Ell. 


Y 


1 noi -j-0fi4 

J-UU J-T^UO^ 


1 


01 


43, 


.3 


05 27 34.8 


n 
u. 


1 fii 


B 


7^0 


900 


Ell. 


Y 


iUiZ-|-UUo 


1 n 


1 o 


20, 


.8 


00 


48 


33 


u. 


. ioO 


r> 


Uo4 


855 


Und. 


Y 


1023-014 


10 


23 


03, 


.9 


-01 24 45.4 


0, 


,150 


B 


2060 


2350 


Und. 


Y 


1047+067 


10 


47 


00, 


.8 


06 45 15.0 


0, 


148 


B 


1176 


1470 


Ell. 


N 


1047-281 


10 


47 


55, 


.3 


-28 07 45.0 


0, 


190 


B 


576 


720 


Ell. 


N 


1151-^117 


11 


51 


15, 


.7 


11 45 10.0 


0, 


176 


B 


366 


456 


Ell. 


Y 


1226+136 


12 


26 


54, 


.6 


13 36 54.0 


0, 


150 


B 


1700 


2250 


Ell. 


N 


1241+095 


12 


41 


10, 


.1 


09 33 31.3 


0. 


,190 


B 


2064 


2580 


Spir. 


N 


1307+085 


13 


07 


16, 


.2 


08 


35 


47 


0. 


,155 


B 


684 


884 


Ell. 


Y 


1325-012 


13 


25 


59, 


.8 


-01 13 47.2 


0. 


.150 


B 


1032 


1290 


Und. 


Y 


1850-782 


18 


50 


08, 


.0 


-78 15 00.0 


0. 


.162 


A 


180 


208 


Ell. 


N 


2140-457 


21 


40 


10, 


.0 


-45 


42 


29 


0. 


,171 


A 


840 


990 


Und. 


Y 



Table 2. The NTT/SUSI2 sample obscrvationnal characteristics, along with host morphologies (spiral, elliptical, or undefined) from 
visual inspection. The presence of signs of interaction is also indicated in the last column. 
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Observation 



Deconvolved Host-QSO 




Figure 3. Left: One of the 3 exposures of the QSO 6QZJ2329- 
2938. Middle; The background and the point source at its fixed 
2 pixels resolution. Right: Residuals of the whole process, i.e. the 
reconvolved model minus the observation. 



ter of each point source, the common numerical background, 
which can be renormalised between different exposures if, for 
example, exposure times differ, and a global shift between 
images. An example of the whole process of simultaneous 
deconvolution is given Fig. [S] It clearly shows that simul- 
taneous deconvolution allows to: (1) seperate efficiently the 
QSO and its host, as the residuals are very good; (2) reveal 
the morphology of the host. 

Among the whole observed sample, only 2 QSOs, namely 
MS13591+0430 and 231711-003604, could not be efhciently 
deconvolved, because they fall very close to the line of sight 
of a very bright star. They were thus removed from the pho- 
tometric analysis, even if visual inspection has given indica- 
tions on their host morphologies. 

We can now use the information provided by the decon- 
volution to compute some physical parameters such as the 
proportion of spiral and elliptical hosts, interacting systems, 
the host and QSO magnitudes in both filters, the center of 
luminosity of the host, its asymmetry coefficient, and try to 
find some correlations between those parameters. 



3.2 Archive data 

As mentioned earlier, 9 QSOs were already observed at high 
resolution, which allow to determine their host morphology 
and magnitudes. If magnitudes are given in the B or R band 
in the original paper, host V m agnitudes are d e duced from 
galaxy colours tables found in iFukugita etHI l| 19951 ). and 
QSO V magnitudes are deduced from apparent V magni- 
tudes taken from the NED database, if available. They are 
upper limits as they include the host luminosity. A summary 
of these additional data is given in Table [1] 



4 ANALYSIS 

Deconvolution of observations obtained in both filters pro- 
vides separated host and QSO images, thus allowing the 
computation of different parameters, which we review one 
by one. 



4.1 Magnitudes 

4.1.1 Apparent magnitudes 

The QSO flux is obtained as an output of the MCS algo- 
rithm. In order to avoid overestimation of the host flux, we 
flagged out the neighbouring objects lying in the field of 
view of the deconvolved image, by visual inspection. The 



host fiux is then computed by summing the intensities of 
the pixels of the deconvolved image. The total flux of each 
component can be converted in apparent magnitudes. This 
requires the knowledge of the apparent magnitude of a stan- 
dard star. Observations of the TPHE B st andard K7-star 
from the Landolt standard stars catalogue l)Landolt|[l993 ) 
allow to infer the zero-point in the V-band, but the appar- 
ent magnitude is not given for that star in the less common 
WB#665 fllter. However, it spans a similar wavelength range 
as the usual R-band. Thus, we deduce the star WB#665 
apparent magnitude by taking the fluxes of Vega and of a 
template K-star spectrum through the R-fllter and through 
the WB#665 fllter: 



rriR — mwB ~ —2.5 log 



f'KBtar*-R 
-F'K3tar*WB 



= -0.03. 



(1) 



Apparent magnitudes in both fllters are listed in Table O 
8 QSOs suffer from significant fiux variations between each 
exposure, indicating non-photometric observing conditions. 
This led us to consider the apparent magnitude taken from 
the brightest exposure in each filter as an upper limit for 
the magnitude (lower brightness limit). 

4.1.2 Absolute magnitudes 

The cosmic recession velocity of the objects we observe im- 
plies a shift in wavelength between the extragalactic emis- 
sion and the observation we make of it. The appropriate cor- 
rection, called K-correction, translate the observed magni- 
tude measured through a filter to the intrinsic emitted mag- 
nitude, given the redshift of the object considered. As the 
host and QSO spectra are different, the conversion of appar- 
ent to absolute magnitudes uses different K-corrections for 
each component. Host galaxies are K-corrected taki ng into 
accou nt their morphologies and redshift according to IPencd 
(| 19761 ). Table 14, which, although rather old, proved to be 
very reliable llKinnev et al.lll996l) QSOs K-corrections are 



taken from ICristiani fc Viol (|l990l ). They are only available 



in the V-band, thus WB#665 absolute magnitudes could 
not be computed. V absolute magnitudes are given in Table 

El 



4.1.3 Errors 

Errors in the computation of magnitudes are dominated by 
two components: the photon noise inherent to the obser- 
vation, and the accuracy of the QSO-host separation. The 
photon noise is estimated by taking the square root of the 
mean QSO+host flux of our sample. That leads to a mean 
error of ~ 0.05%. 

To estimate the error on QSO-host separation, we first 
selected 3 representative QSOs with different quality de- 
convolution residuals: 1241-1-095, HE0146-3755, and Q2252- 
2434. For each of these QSOs, we varied the value of the 
QSO intensity until clear signs of bad PSF subtraction ap- 
pear (clear hole or peak at the QSO position). Signs of bad 
PSF subtraction begin to appear when we vary the QSO 
intensity by 1.1 to 7.0 per cent, suggesting that this is the 
level of uncertainty in the QSO magnitudes. Propagating 
these errors to magnitude errors leads to aM — 0.02 — 0.05. 

Another way to test the error bars on galactic and QSO 
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Magnitudes Asymmetry 



Name 




QSO 








Host 




"■QSO 






my 


M 


V 




my 


My 




V 


WB 


V 


WB 


005709+144610.1 


15.31 


-24 




15.26 


16. 


,95 


-22.49 


16.34 


n sA 


n Q9 
u.yz 


n 9R 
u.zu 


u.o 


011110-101631.8 


15.94 


-24, 


.02 


15.81 


16. 


.87 


-22.84 


16.29 


0.7 


0.64 


0.27 


0.21 


011845+133327.1 


17.49 


-22 


.oo 


17.24 


16, 


,97 


-22.97 


16.28 


n 1 


12 


14 


n 1 


015530-085704.0 


16.26 


-23 


.46 


16.16 


17, 


,02 


-22.6 


16.55 


0.23 


0.91 


0.38 


0.53 


021218-073719.8 


17.16 


-22 


. 1 u 


17.17 


17, 


,55 


-21.84 


17.19 


u.uo 




12 


12 


021360+004226.7 


14.55 






- 


19, 


,54 


- 


- 


1 '^A 




n '^A 


n 1 Q 


025007+002525.3 


16.93 


-zo 


44 


16.87 


18, 


,69 


-21.57 


18.01 


n Qi 


Ql 
u.y J- 


n 9"^ 
u.zo 


VJ. lO 


032214+005513.4 


15.6 


-24 


.71 


15.52 


17, 


,87 


-21.92 


17.27 


0.64 


0.77 


0.15 


0.35 


101044+004331.3 


15.6 


-24 


.oo 


16. 


17, 


,52 


-22.17 


16.7 


u.uu 


n 1 R 

U. J-U 


U. AU 


n (17 

U.U i 


113706+013947.9 


16.22 


-23, 


.89 


16.23 


18. 


,34 


-21.52 


17.64 


0.15 


0.65 


0.26 


0.8 


122534-024757.2 


16.93 






16.85 


18. 


,39 


- 


17.86 


41 


u.uy 


n 1 R 
u. ±u 


22 


161532-002730.3 


16.7 


-23, 


.00 


16.47 


16, 


,71 


-22.57 


16.43 


0.1 


0.53 


0.08 


0.08 


205032-070131.2 


16.95 






16.96 


16, 


,75 


- 


16.1 


24 


u.oy 


U. J-O 


u.zo 


232260-005359.3 


16.73 






- 


16, 


,56 


- 


- 


0.74 


0.71 


0.45 


0.17 


235156-010913.3 


14.55 






- 


19, 


,54 


- 


- 


1 .oo 


n 89 
u.oz 


1 25 




Q 0022-2044 


15.59 


-24 


.19 


15.18 


18, 


,59 


-20.73 


20.52 


0.73 


1.45 


0.23 


0.33 


CT 289 


16.41 


-23 


.21 


16.38 


18, 


,71 


-20.47 


17.97 


0.24 


0.24 


0.13 


0.22 


MS 01325-4151 


17.54 


-22 


.25 


17.21 


17, 


,53 


-22.16 


16.94 


0.30 


0.43 


0.38 


0.22 


MS 10302-2757 


15.35 


-24, 


.24 


15.40 


16. 


.79 


-22.38 


16.36 


0.47 


0.68 


0.23 


0.05 


PKS 1241-399 


17.96 


-22, 


.43 


17.74 


17. 


.67 


-22.60 


16.88 


1.02 


1.1 


0.65 


0.64 


Q 1421-0013 


15.65 


-23 


.94 


15.64 


17, 


,47 


-21.91 


17.64 


0.54 


0.29 


0.13 


0.23 


PDS 456 


13.42 


-28 


.20 


13.10 


16, 


,91 


-24.60 


16.30 


1.66 


2.3 


1.4 


0.36 


MS 20078-3622 


16.27 


-zo 


.ol 


16.30 


17, 


,73 


-21.85 


17.48 


42 


u.uo 


7"^ 
U. i o 


n 71 

U. 1 -L 


Q 2240-2411 


17.49 


-22 


.50 


17.39 


16, 


,73 


-23.02 


16.22 


0.87 


0.84 


0.18 


0.28 


Q 2252-2434 


15.44 


-24 


.01 


15.55 


16, 


,89 


-22.47 


16.15 


0.97 


0.87 


0.33 


0.14 


6QZ J232927-2938 


15.85 


-24 


.25 


15.82 


17, 


,83 


-22.00 


16.98 


0.43 


0.21 


0.14 


0.05 


HE 0027-3118 


16.10 


9*^ 


9Q 
.zy 


16.19 


17. 


,53 


-21.44 


16.88 


n ^9 


42 


12 


n 9fi 
u.zu 


HE 0108-5422 


16.70 


-zo. 


.O J. 


16.62 


19. 


.01 


-20.47 


17.65 




u.oo 


n 9Q 

u.zy 


U.U'l 


HE 0146-3755 


15.98 


-23 


.43 


16.03 


16. 


.38 


-22.83 


15.89 


0.72 


0.47 


0.21 


0.35 


HE 0226-5209 


16.64 






17.26 


16, 


,84 


- 


16.74 


1 12 


1 21 


n RR 
u.uo 


U.Oo 


HE 0227-4123 


18.84 






- 


19, 


,59 


- 


- 


0.24 




.37 




HE 0250-4400 


15.57 


-24 




15.56 


17, 


,35 


-21.92 


17.06 


n 7 




21 


n 07 

U.U / 


HE 0441-2826 


14.45 


-25 


.16 


14.32 


16, 


,04 


-23.13 


16.29 


0.73 


1.09 


0.08 


0.11 


HE 1101-0959 


17.08 


-23 


.03 


17.04 


17, 


,69 


-21.89 


17.26 


0.39 


0.47 


0.08 


0.08 


HE 1202-0501 


16.15 


9*^ 


.U ( 


16.06 


16 


;.5 


-23.21 


16.28 


1 62 


1 7*^ 


1 23 


1 22 


HE 1211-1905 


16.45 


9*^ 




16.55 


16, 


,64 


-22.70 


16.13 


1 14 


1 4 


41 


RO 
U.UU 


HE 1236-2001 


16.04 


-24, 


.22 


15.85 


17. 


.21 


-22.50 


16.75 


0.51 


0.46 


0.21 


0.16 


HE 1255-0437 


16.34 


9*^ 


.^y 


16.35 


16. 


,86 


-22.50 


16.25 


u. ±o 


n 9Q 

u.zy 


n "^9 
u.oz 





HE 1256-2139 


15.79 


-23, 


.93 


15.75 


18. 


,81 


-20.49 


18.75 


1.9 


1.68 


0.43 


0.61 


HE 1300-0657 


16.53 


9*^ 




16.48 


16. 


.82 


-22.94 


16.22 


n ^9 
u.oz 


n 97 

U.Z 1 


u. ±o 


04 


HE 2345-3939 


16.61 


-23 


.50 


16.65 


16, 


,37 


-23.47 


15.83 


1.17 


1.14 


0.47 


0.38 


0056-363 


14.83 


-24 


.ol 


14.9 


17, 


,31 


-21.9 


19.27 


1 49 


J- .oo 


1 

1 .oo 


_L .O^ 


0132+077 


17.37 


-22 


1 S 


17.4 


16, 


,96 


-22.39 


16.32 


n 9Q 
u.zy 


n 9Q 

u.zy 


n OR 
u.uu 


OS 
u.uo 


0213-484 


16.92 


99 

-zz. 




16.92 


17, 


,67 


-21.64 


17.09 


U.OD 


n 1 7 




n 1 o 
U.iy 


0357+107 


15.81 






-24.88 


18. 


,21 


_ 


-21.98 


0.85 


0.71 


0.06 




1001+054 


15.57 


-24, 


.12 


15.69 


17. 


,76 


-21.85 


17.31 


0.53 


0.8 


0.24 


0.37 


1012+008 


15.15 


-24, 


.89 


15.15 


16. 


,45 


-23.47 


15.76 


2.08 


1.95 


0.85 


0.6 


1023-014 


18.85 


-20, 


.73 


18.89 


17. 


,49 


-22 


17.38 


2.04 


1.95 


0.51 


0.18 


1047+067 


16.82 


-22, 


.65 


16.84 


17. 


.12 


-21.82 


16.49 


0.09 


0.34 


0.09 


0.12 


1047-281 


15.27 


-24, 


.99 


15.31 


17, 


,54 


-22.18 


16.76 


0.28 


0.33 


0.1 


0.48 


1151+117 


15.61 


-24, 


.28 


15.64 


18, 


,79 


-21.00 


17.59 


1.51 


1.37 


0.4 


0.32 


1226+136 


17.43 


-22 


.08 


17.07 


17, 


,06 


-22.02 


16.61 


0.18 


0.38 


0.26 


0.22 


1241+095 


16.81 


-23, 


.24 


16.79 


17, 


,06 


-22.74 


16.73 


0.74 


0.65 


0.18 


0.07 


1307+085 


14.79 


-24, 


.84 


15.29 


17, 


,42 


-21.78 


17.75 


1.29 


0.28 


0.29 


0.15 


1325-012 


16.95 






16.96 


17, 


,24 




16.36 


1.59 


1.55 


0.19 


0.18 


1850-782 


14.99 


-25, 


.20 


14.94 


17, 


,00 


-22.75 


16.35 


0.42 


0.55 


0.13 


0.1 


2140-457 


15.68 


-24 


.12 


15.64 


17, 


,92 


-21.77 


17.47 


0.45 


0.39 


0.12 


0.12 



Table 3. Apparent and absolute magnitudes arc given separately for the QSO and the host in the V-band, while only apparent 
magnitudes could be computed in the WB#665-filter. For non piiotometric observations, absolute V-magnitudes are not displayed and 
apparent magnitudes must be taken as a lower limit. Asymmetries computed with respect to the QSO center and to the center of 
luminosity of the host are given in both filters (see text for details) . 
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Regular 


Interaction 


Multiple nuclei 


Total 


Spiral 


17 


8 





25 


Elliptical 


23 


28 


4 


51 


Undef. 





24 


4 


24 


Total 


40 


60 


9 





Table 4. Occurcncy (in % of the whole sample) of morphology 
classes and interaction features. 



magnitudes obtained through deconvolution is to add a sim- 
ulated QSO on an observed galaxy in one of the observed 
fields and run the deconvolution process on this new image. 
The input and output fluxes may then be compared to esti- 
mate the error bars on magnitudes. After running this pro- 
cedure for several objects, we observe a deviation of at most 
20% for the host flux and 2% for the QSO flux. Not surpris- 
ingly, the uncertainty for the host is stronger when the QSO 
intensity is high relative to the host. On average, given the 
nuclear/host ratios in the sample, we can consider an un- 
certainty of 0.1 mag for the host magnitude and 0.01 for 
the QSO magnitude. The errors coming from the QSO/host 
separation thus clearly tend to dominate the photon noise. 



4.2 Morphology 

One of the first natural steps in the study of QSO-host sam- 
ples is to seek if there is a link between the host morpho- 
logical type and the QSO activity. More precisely, the cur- 
rent open questions are: Do gravitational interactions trig- 
ger QSO activity? If yes, what are the physical processes 
at work? Are there other ways of triggering activity? Which 
ones? Once the nucleus is turned on, how does it evolve with 
its host?... One way to answer these questions is to find clas- 
sification criteria for host galaxies that enable to compare 
them with quiescent galaxies. The bimodal distr i butio n of 
galaxy colour-magnitude diagrams tealdrv et all \2004) al- 
lows such a comparison. Namely, 'Marti n et al.l (|2007l ) find 
that the host galaxies of AGNs lie more often in between the 
red (corresponding to early-types) and blue (corresponding 
to late-types) sequence of such diagrams, suggesting that 
the AGN phenomenon is a transitionary step between those 
sequences. Another possibility is to simply classify hosts by 
their morphological type and find out the frequencies of in- 
teracting, elliptical and spiral systems, in which QSOs are 
found to lie, and compare them with previous results. 
We chose to classify each host galaxy according to its mor- 
phology by visual inspection. Namely, the different selection 
properties are: (1) Spiral, elliptical or undefined morphol- 
ogy (exclusive), where the undefined morphology most often 
corresponds to cases in which the merging process is violent 
enough to totally disturb the host, preventing the deter- 
mination of its morphology, and, (2) signs of gravitational 
interaction and presence of more than one point source (non- 
exclusive). The percentages, calculated on the total sample 
(SUSI2-I- archive sample), can be found in Table ID 

The proportion of ellipticals (51%) and spirals (25%) 
is essentially comp ati ble with some previo us results, such as 
IShade et all ()2000f ) or lLetawe eraD (|2007[ ). Ho wever, we find 



a lowe r proportion of elfiptical gala xies than in Dunlop et al.l 
' 51), (Hamilton et~al] (|2008l) and iFlovd et al.l l|2004l) . but 



this may be due to a se l ection effect in their sample, as men- 
tioned in lLetawe et al.l l|2007l ). Concerning the proportion of 
hosts showing signs of gravitational interaction, we just men- 
tion for the moment that the 60% we find are much higher 
than the 5 to 20% found for quiescent galaxies in some pre- 
vious studies dLe Fevrd 120001: iBell et al.ll2005l : [Zheng et al.l 
2004), although IShi et al.l l|2009l ^ argue in favour of a higher 
merger rate. One difficulty in estimating the proportion 
of merger remnants is that, as shown by simulations, the 
timescales during which they show clear signs of disturbance 
are often shorter than the total merger timescale (|Lotz et al.l 
2008). On the othe r hand, we find a proportion of interac- 
tion very similar to lSchmittI (|200ll) for both active and inac- 
tive galaxies. Namely, ~ 55% of all our elliptical hosts show 
signs of gravitational interactions, and this proportion drops 
to ~ 32% for spirals, while they find respectively ~ 50% and 
~ 25%. If we consider our undefined class as strong merg- 
ers, o ur percentage of 24% is very similar to I Greene et al.l 
( |2009l ). who find a proportion of one quarter for highly dis- 
turbed morphologies in obscured active galaxies. 

In conclusion, we can state that our estimates show that 
the morphologies of active galaxies are not different from the 
quiescent ones, and that the occurence of merger signs in 
QSO hosts is higher or similar to inactive galaxies, regarding 
previous discordant studies. If higher, we could assess the 
infiuence of mergers on activity. Thus detailed studies on 
the presence of merger signs in inactive galaxies at similar 
spatial resolution would be required before one can firmly 
conclude. 

All in all, the major difficulty in comparing different 
studies lies in the definition and quantification of the degree 
of interaction in a system. It is often done by simple eye- 
check, which is always subjective. That leads us to define 
the following asymmetry coefficient. 



4.3 Asymmetry 

An idea that has been recurrent for the past 10 years is to 
quantify interactions and mergers via the asymmetry of the 
system. Several methods can be found in the literature. We 
briefiy review the most releva n t ones and their results. First 
of a ll. lAbraham et al.l ljl994f ): IConselice et"al] (|2000l . 120031 ) 
and IShi et al.l l|2009D used an asymmetry index defined es- 
sentially as the subtraction of the image of a galaxy by the 
same image rotated by 180°. This simple method proved 
to be a good tracer of the degree of interaction in merging 
galaxies. With the help of two supplementary parameters, 
t he co ncentration index and the dumpiness. Consel ice et al.l 
(120031 ) achieve to describe all usual morphological types in a 
totally quantitive way, solely from image analysis. Another 
method called Zurich Estimator of Structur al Types (ZEST) 
has been proposed bv lScarlata et all (|20o3). It uses the Ser- 
sic index n returned from a galaxy fit and 5 other basic 
nonparametric diagnostics (amongst which the asymmetry 
coefficient) to quantify the properties of galaxy structure, 
using Principal Component Analysis. They find strong evo- 
lutionary effects between z — and z = 0.7 for faint galaxies 
{Mb > -20.5). 

Another morphological parameter used as a signpost of a 
nonequilibrium global dynamical state is the so-called lop- 
sidedness. It is defined as the radially averaged m = 1 az- 
imuthal Fourier amplitude measured between radii enclos- 
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ing 50% and 90% of the galaxy light. Using this definition 
for a low-redsh i ft sam ple (z < 0.06) of ^ 25000 galaxies, 
iReichard et all (|2009l l find a trend for more powerful AGN 
to be hosted by more lopsided galaxies, and conclude that 
is due to the strong link of the age of the stellar population 
in the galaxy bulge to both the AGN luminosity and to the 
lopsidedness of the host. 

These analyses concentrate on Type 2 AGNs or quiescent 
galaxies, and do not take into account Type 1 QSOs, which 
considerably overweight the central regions of the host in 
the computation of the asymmetr y index. In order t o in- 
clude QSOs in that kind of analysis, iGabor et all (|2009h first 
separe QSO and host components with the GALFIT tool, 
and then compute the asymmetry index in the same way as 
IConselice et all l|2000l l. They find that, at 0.3 < z < 1.0, 
QSO host asymmetries are no more prevalent than in qui- 
escent galaxies. Another way of co mputing as y mmet ries for 
QSO host galaxies is proposed by iKim et al.l (|2008l ). using 
an all-GALFIT computation. After separing the QSO and 
the host, Fourier components are fit to the images showing 
significant nonaxisymmetric features. Their result is that, 
for low redshift QSOs (z < 0.35), galaxy mergers and tidal 
interactions seem to play an important role in regulating 
and fueling the nuclear activity. The apparently opposite 
resu lts in different re dshift ranges from iGabor et al.l (|2009l ) 
and iKim et aP (|2008l ) are suggestive of a cosmological evo- 
lution of the importance of interactions. It reinforces our 
belief that the study of asymmetries is a powerful tool to 
understand QSO-host interactions. 

Here, in order to get rid of any assumption concerning the 
host morphology or the estimated center of the host, we pro- 
pose an alternative way to define the asymmetry, inspired 
by the statistical third order moment, called skewness in- 
dice. More precisely, if I{x) is the intensity in the pixel x, 
we define the center of luminosity of the host by 



(2) 



where the sums are over the total number A*' of pixels in a 
square subimage defined by the user and containing the host 
galaxy. 

Then, we define an asymmetry coefficient with respect to a 
center c located at some position Xc 



(3) 



This method has the advantage of being totally model 
independent, and takes into account the whole light from 
the host. It thus requires no a priori knowledge about the 
shape of a galaxy, which involves the risk to misguide any 
interpretation. 

Another advantage of this definition is that c can be chosen 
to be either the center of luminosity cl, thus describing the 
asymmetry of the galaxy alone, or the QSO position cqso, 
thus giving information about the asymmetry of the whole 
system (QSO-f-host). From now on, unless clearly stated 
otherwise, the term "asymmetry coefficient" will refer 
to the asymmetry with respect to cqso- Indeed, a QSO 
located off-center in an otherwise seemingly symmetrical 
galaxy can be considered as a sign that something special is 
happening. For example, in the framework of binary black 
hole mergers, the resulting black hole is supposed to recoil 
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Figure 4. Asymmetry coefficients in the V-band a(V) are dis- 
played for 6 QSO hosts which are at different levels of interaction. 
A higher asymmetry value corresponds to a higher degree of in- 
teraction. The spatial scale is identical in all frames. 



in a direction opposite t o the gravitational wave emission 
(jRedmount fc Ree3ll989l ). During this recoil, the accretion 
disc may stay bound to the bl ack hole and consequent l; 
shine off-centered. Simulations ( Volo nteri fc Madaul 1200. 
predict that a population of off-nuclcar AGNs may already 
be detectable at low and intermediate redshifts. 

In order to minimize the importance of border effects 
linked to the choice of the square subimage where the 
asymetry coefficient is computed, we multiply the region of 
interest by a broad Gaussian function, centered on the center 
of luminosity. A FWHM of 30 pixels, matching the extent 
of the galaxy, reveals to be appropriate in the majority of 
cases as it is not too narrow, which would give more weight 
to asymmetries due to shifts between and cqso in com- 
parison with galactic asymmetries, and not too broad to be 
useless because including unrelated features in the field. The 
FWHM was only changed when the galaxy was significantly 
smaller or larger (the extreme values for the sample being 
down to 15 and up to 45 pxl). If the width of the Gaussian 
is suited to the galactic extend, a difference of 2pxl FWHM 
does not induce significant changes in the asymmetry. In 
Fig. 131 six deconvolved images of QSOs and their hosts are 
shown along with their asymmetry coefficient. This illus- 
trates the relevance of using the asymmetry coefficient as a 
measure of the degree of interaction of a system. 

The global sample can be divided into two subsambles, 
one with a{V) ^ 0.7 (low asymmetry), and the other with 
a{V) > 0.7 (high asymmetry). This value seems appropriate 
(see Fig. for a good separation between minor interaction 
events and galactic-scale mergers involving several galaxies. 
The low-asymmetry subsample contains 60% of the sample. 
Let us mention that asymmetry coefficients were not com- 
puted for the archive data because hosts were not separated 
from their nucleus. They are thus not included in forthcom- 
ing analysis including asymmetries. 



Errors 

Errors on the asymmetry coefficient and cl have been esti- 
mated by varying the input parameters (size of the square 
subimage in which we compute the coefficient and Gaussian 
FWHM) and see how it influences the result. First, it shows 
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that cl is very stable (a 10% variation of the subimage size 
leads to a < 1% variation of cl). Conversely, the asymmetry 
coefficient is quite sensitive to the Gaussian FWHM (a 25% 
variation in FWHM leads to a ~ 15% variation in asym- 
metry coefficient). However, the restrictions for setting the 
FWHM value explained above do not allow such large vari- 
ations, and thus we do not expect a change of more than 
> 10% in the asymmetry coefficient, which is clearly suffi- 
cient for the present purpose. 

Two other sources of error must be taken into account: the 
variation in sampling due to the redshift range spanned, 
and S/N variations, which might change the influence of the 
background noise on the asymmetry. 

First of all, the exposure times were chosen to reach a com- 
parable S/N for all the hosts, thus considerably lowering the 
risk of error due to S/N variations. Secondly, it might be ex- 
pected that higher-z hosts would show a lower asymmetry 
than their lower-z counterparts, as the hosts substructures 
are less and less resolved with increasing redshift. In order to 
test the importance of this effect in our sample, we use the 
PYRAF tool "magnify" to decrease the resolution of the ob- 
servations of the most asymmetric hosts {a{V) ~ 0.6) among 
the nearest objects {z < 0.15) from around 0.4 kpc/pxl to a 
0.52 kpc/pxl scale, corresponding to the upper limit of the 
sample, z = 0.198, and compute the asymmetry on this new 
image. The result is that the asymmetry coefficient drops 
at most by 0.02, revealing that the redshift range is small 
enough to avoid such biases in the estimates. We can thus 
neglect both of these effects. The calculated asymmetries 
with respect to cqso and c_l are given in Table O for both 
fllters. 



5 RESULTS 

Now that the parameters have been introduced, we aim to 
flnd correlations between them and try to specify some as- 
pects of QSO-host interactions. 
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Figure 5. My (Host) is plotted versus My{QSO). Solid line: 
the b est fit relation of our sample. Dashed line: fit from 
[Hamilton et al.l (|200_8). Dotted line: fit from Dunlop ct al. (2003|). 
Filled symbols represent non-interacting galaxies, while open ones 
are for galaxies showing signs of interaction. Triangles are for spi- 
rals, hexagones for ellipticals, whereas squares stand for unclassi- 
fied morphologies. 



Subclass 


a 


b 


P 


All 


0.21 ±0.10 


-22.20 ±0.11 


0.28 


Spirals 


-0.28 ±0.39 


-22.35 ±0.26 


0.20 


Spirals w/o inter. 


-0.5 ±0.61 


-22.26 ±0.32 


0.26 


Ellipticals 


0.23 ±0.14 


-21.94 ±0.13 


0.31 


Ellipticals w/o inter. 


0.17 ±0.17 


-22.02 ±0.20 


0.29 


Interactions 


0.26 ±0.12 


-22.22 ±0.14 


0.39 



Table 5. Linear Mv(Host) = a{Mv(QSO) - {MviQSO))) + b 
relation for different morphological subclasses, p is the Pearson's 
correlation coefficient. 



5.1 QSO-Host absolute magnitudes relation 

A common question in QSO hosts studies is to seek 
if there is a correlation between QSO and host abso- 
lute magnitudes. Theoretically, such a trend is expected 
as a consequence of the more fundamental and well 
constrained M ass(Spheroid)- Mass(Black Hole) r elation 
ijFerrarese fc M orritt 2000; Mar coni and Hundlioosh . How- 
ever, it might be strongly disturbed by the variety of accre- 
tion rates found in AGN, which prevents from linking the 
black hole masses and magnitudes unequivocally. Equally, 
spheroid masses are not expected to match the galaxies 
magnitudes, as the latter contain significant supplementary 
structu res (tidal tails, spiral arms and bars, discs). Previous 
resul ts (|Hamilton et al.lboOStlFlovd et al.ll2004l:IShade et al.l 
|200(|) tend to indicate only a slight correla tion, while oth- 
ers teahcall et al.l[l997l : IPunlop et all 120031 ') find a relation 
compatible with no correlation at all. Our sample shows a 
weak correlation between host and QSO magnitudes (Fig. 
O where our correlation is compared with mean relations 
from other samples, and the morphology classification is in- 
dicated). Namely, we find that, if < Mv{QSO) > is the 



mean value of the whole sample, 

Mv{Host) = (0.21 ± 0. 10) (Mv(Q5'0)- < My {QSO) >) 
- (22.20 ±0.11), (4) 

with a Pearson's correlation coefficient p = 0.28 (the corre- 
lation is calculated with MviQSO) - {My (QSO)) instead 
of MviQSO) so that the intersection with the Mv(Host)- 
axis corresponds to the average case where MviQSO) = 
(MviQSO))). 

In Table O the relations found for different subclasses 
are presented with associated correlation coefficients, all 
globally weak. Brighter ellipticals most probably harbour 
brighter QSOs, as theoretically expected, but this does not 
hold for spirals, resulting in a poor correlation between 
MviHost) and MviQSO) for the whole sample, in agree- 
ment with the aforementioned previous studies. 

As explained above, the asymmetry coefficient is a good 
indicator of the degree of disturbance due to interaction in 
a system. Thus, we can replace the traditional morphology 
classification used in Fig. [5] by an asymmetry-based classifi- 
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Subclass 




a{y) 




(j{WB) 


Spirals 


0.59 


0.36 


0.58 


0.39 


Ellipticals 


0.57 


0.45 


0.61 


0.42 


Interactions 


1.11 


0.66 


1.23 


0.63 



Table 6. Average asymmetries are given in the two filters for 
spirals, ellipticals, and hosts with signs of interaction. 
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Figure 6. Top. Mv{Host) is plotted versus MviQSO) with the 
QSOs labelled according to their asymmetry subsample. Bottom. 
Deviations AM from the best fit are plotted against the asymme- 
try coefficient, with a linear fit on the postive and negative values 
of AAf separately. 
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cation. Figure [S] shows the My {Host) — MviQSO) relation, 
with this new subsample separation. 

We find that the high-asymmetry subsample shows a 
dispersion twice larger around the fitted relation. This trend 
is even clearer when we plot the deviations from the fit as a 
function of the asymmetry (Fig. |6l bottom) . It shows a cor- 
relation between both variables indicating that the higher 
asymmetry hosts deviate more strongly from the fit. 
Moreover, the asymmetric hosts tend to lie on average below 
the fit in Fig. |6l It means that, for a given nucleus magni- 
tude, highly asymmetric systems tend to have brighter host 
galaxies. It can simply be explained by the fact that impor- 
tant mergers, corresponding to high a{V) values, contain 
two or more galaxies which boost the value of the computed 
magnitudes. More interesting are the systems which have a 
surprisingly low-luminosity host (or powerful nucleus) along 
with highly asymmetric features (open squares at the top of 
Fig- Ell - Most probably, they are systems where interactions 
have most strongly enhanced th e QSO activity and l uminos- 
ity, as suggested for example in iLetawe et al.l (|2007l) . 

It is also interesting to check if asymmetric hosts in a 
filter are also asymmetric in the other one. Figure [7] shows 
that for the vast majority of the QSOs observed, asymme- 
tries are, as expected, comparable in both filters. The best 
fit is a{WB) = (0.86 ± 0.08)a(V) + (0.15 ± 0.07), with a 
correlation coefficient p = 0.85. However, some objects lie 
way out of the mean relation. They will be discussed in- 
dividually in the Appendix. Figure [7] is useful to describe 
the link between the two classification methods (morphol- 
ogy vs asymmetry). It is clear that the most disturbed cases, 
morphologically undefined and indicated by an open square, 
also have the highest asymmetry value. Spirals and ellipti- 
cals have comparable asymmetries, whatever the filter. Ta- 



Figure 7. Correlation between the asymmetry coefficients mea- 
sured in both filters. Note that all the high asymmetry sys- 
tems have also been classified as interacting systems. Filled sym- 
bols represent non-interacting galaxies, while open ones are for 
galaxies showing signs of interaction. Triangles are for spirals, 
hexagones represent ellipticals, whereas squares stand for unclas- 
sified morphologies. 

ble |6] summarizes the mean asymmetries in both filters as a 
function of morphology. 

5.2 Colours 

As explained in Section 2, one of the goals of the study is 
to compare the stellar and gas distributions in order to seek 
for cases where there are significant differences between both 
components. We thus take the difference in apparent mag- 
nitudes between both filters, the colour Am — my — mwB, 
for the host and the QSO separately. As a preliminary com- 
ment, let us mention that the deconvolution process, as it 
treats both filters exactly in the same way for separating the 
QSO and host components, should have no influence on the 
derived colours. Furthermore, during the deconvolution pro- 
cess, given the slight uncertainty on magnitudes induced by 
QSO-host separation as mentionned in Section 14.1.31 some 
flux from the QSO might be interpreted as galaxy flux. The 
QSO being bluer than galaxies, it would lead to an overes- 
timation of hosts blue colour compared to normal galaxies. 
Such a perturbation is however hard to estimate. It might 
lead to a slight shift of our hosts distribution in the his- 
togram towards the redder galaxies in both filter. It will 
nevertheless not change the conclusions obtained from the 
colour analysis. 
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Figure 8. Histogram of Am{Host) with 0.2 mag bins. Typical values for template galaxies from starburst to E types are indicated 
by the shaded strip. Hosts with the strongest emission lines appear to the left. Left: global sample. Middle: low and high-asymmetry 
subsamples, which have a different shape, the high asymmetry hosts being shifted to bluer colours, meaning that they display both 
stronger emission lines and/or a bluer continuum. Right: ellipticals, spirals, and undefined morphology galaxies. Elliptical hosts have the 
broadest range of colours. 



First of all, one might expect QSOs to be brighter 
in the V-filter, as it contains the strong H/3 and [OIII] 
emission lines. We find a mean colour for the QSOs of 
Am{QSO) = 0.03 ± 0.17. In order to check if this is 
not due to any proces sing error, the spe ctrum of a typi- 
cal QSO HE1302-1017 l|Letawe et al.ll2007l ) {z = 0.278) has 
been blueshifted to z = 0.16 and integrated in both filters. 
Converting the flux ratios into magnitude differences leads 
Am(HE1302-1017) = 0.08, which lies well within our error 
bars. 

For the host galaxies, we find an average Am(host) — 
0.41±0.58. The left panel of Fig. [8] shows the distribution 
of Am for all our host galaxies. Some of these host galax- 
ies have very blue colours Am < 0, even bluer than the 
QSOs. To identify the origins of the blue light, we compare 
the colours of the sample with the colours of galaxy tem- 
plate spectra, estimated as for QSO colours in integrating 
spect ra in both filters. We used galax ies of types E to Sc 
from iMannucci et al.1 (|200ll) . and from lKinnev et all (|l996l ) 
which constructed 6 starburst galaxies with different extinc- 
tion values (from SBl to SB6 with increasing extinction). 
The associated colours are reported in Table [7] together with 
the colours for our sample. 

The colour is directly influenced by the strength of the 
ionization lines, but also by the slope of the continuum. 
First, we estimate an upper limit on the colour variation 
due to the continuum contribution by removing the ion- 
ization lines from the templates and compute Am again. 
The values we obtained are presented in the third column 
of Table [7] It shows that the colour associated to the bluest 
slope for the continuum is estimated at 0.34. Moreover the 
slope of the continuum may be responsible for a variation 
in Am which should not go beyond 0.5, which is the differ- 
ence between the reddest (E type) and bluest (SBl type) 
spectra. On the other hand, placing upper limits on the 
emission lines contribution is harder, as QSO hosts might 
have stronger lines than the templates spectra. The host of 
HE1434-1600 is a good example of this as it has particu- 
larly strong emission lines wit h a relatively weak and red 
continuum (|Letawe et al.|[20o3 ). From the spectrum of this 



Galaxy Am Am Diff. Am 

Global Continuum Sample 

E 0.87 0.87 Ellipticals 

SO 0.86 0.86 0.27±0.78 

Sa 0.85 0.85 

Sb 0.76 0.77 -0.1 Spirals 

Sc 0.65 0.67 -0.02 0.60±0.15 

SB6 0.5 0.54 -0.04 

SBl 0.23 0.34 -0.11 Undefined 

HE1434-1600 0.27 0.63 -0.36 0.58±0.28 



Table 7. Colours estimated on galaxy template spectra and on 
the peculiar QSO host H E1434-1600 with a low, red continuum 
and strong emission lines llLetawe et al.ll2004h . Colours are given 
with and without the contribution of emission lines, together with 
the difference in colours between both (total minus continuum 
only). The last column gives the average values and scatter for 
our QSO host galaxies as a function of morphology 



QSO, we compute Am{Total) = 0.27 and, after removal of 
the ionization lines, Am{C ontinuum) = 0.63. So, the ion- 
ization lines may induce variations in Am of at least 0.36. 
In conclusion, both contributions influence the result at a 
relatively similar level. For colours lower thant 0.34, ionized 
gas must be present, and for variations from templates larger 
than ^ 0.5 a bluer than expected colour cannot be explained 
only by young stellar population, also requiring the presence 
of ionized gas. 

In order to analyze the influence of morphology and 
asymmetry on the colour distribution, we show on Fig.|8]his- 
tograms of the colours for different hosts subsamples (global 
sample and subsamples with different asymmetry or mor- 
phology, from left to right). We plot the numbers of hosts 
whose Am is contained in bins of 0.2 mag, from —2. to 1.5, 
and compare them to the template galaxy colours mentioned 
above. 

First, Fig. [8] shows that the range spanned in Am for 
our hosts is much larger that the range spanned by the tem- 
plates, only covering the 0.3 — 0.9 region. Moreover, the ma- 
jority of the hosts seems to have colour typical of Starburst 
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to Sc galaxies, characterized both by a large amount of ion- 
ized gas and a blue stellar continuum. 

Secondly, there is a tendency for the most asymmet- 
ric hosts to be bluer. Even if the dispersion is quite large 
for the high asymmetry sample, (Am(high-a(V) hosts) — 
0.17 ± 0.81, whereas Am(low-a(V) hosts) = 0.61 ± 0.24, 
leading to uncertainties on the mean values of resp. ±0.1 and 
±0.04 given the size of these subsamples), we notice that the 
bluest host are highly asymmetric. These objects are inter- 
preted as hosts involved in major mergings, having a higher 
proportion of ionized gas and/or a bluer continuum than 
non-interacting ones. A sim ilar tendancy has been found in 
iLetawe et all l|2007l [ioOSbh . where links between gas ionized 
by the QSO radiation and gravitational interactions are ob- 
served. Thus, it is tempting to propose that our observations 
give a statistically more relevant proof of this kind of inter- 
relation between QSOs and their hosts. However we cannot 
disentangle the different ionization processes, namely ioniza- 
tion by the QSO, by shocks, or by stars. Consequently, the 
higher degree of ionization in highly asymmetric hosts might 
also be interpreted as due to prominent HII regions typical 
of strong star formation, enhanced during the interaction 
process. 

Third, we point out the fact that all hosts that have 
Am{Host) < 0, thus the bluest ones, are ellipticals. It is 
a suprising result as ellipticals are often supposed to be 
evolved galaxies with old stellar populations and only a 
low proportion of gas and dust. This trend for elliptical 
QSO host galaxi es to have bl uer colours was already re- 
ported by Jahnk e et al.l I 2004al) for similar ob jects, and by 
ISanchez et al.l (|2004l ) or 



Jahnke et al.l l|2004bl ) for samples 



at higher redshift, and was explained by a younger stel- 
lar component in the h ost. Spectral analysi s of a sample 
of nearby QSO hosts in iLetawe et al.l ([20o3) also pointed 
to a younger-than-average stellar population. As the colour 
variation due to the continuum (thus the stellar popula- 
tion) presented above reaches at most 0.5 mag, and as the 
bluest slope only leads to Am=0.34, emission line contribu- 
tion has to be present in at least a dozen of outliers having 
Am{Host) < 0. They likely correspond to a class similar to 
the alr eady studied QSO HE1434-1600 (|Letawe et al.ll2004l . 
l2008b[ ). which has a seemingly normal elliptical host, but 
also filamentary structures consisting of gas ionized by the 
central nucleus on both sides of it. This system is in grav- 
itational interaction with a smaller neighbouring elliptical. 
Indeed, if the ionization was due to stars in those elliptical 
outliers, it would be expected to be at least as strong in spi- 
rals, which is not the case. Thus, if an ionization source has 
to be favoured, QSO or shock ionization clearly seems to be 
preferred. 

We finally see in Fig. |8] (right panel) that spirals 
better match with classical values (inactive galaxies) than 
ellipticals, the latter being responsible for the broadening 
of the range spanned by the whole sample. 

This hosts colour analysis reveals the difference between 
QSO hosts and quiescent galaxies. QSO hosts tend to have a 
bluer continuum, but also to contain more ionized gas than 
their quiescent counterparts, and this trend is stronger in 
ellipticals and highly asym metric systern s . It i s consistent 
with the previous study of lScoville etall l|2003l ). who find, 
by analysing CO emission of 12 low redshift {z < 0.1) PC 



QSOs, that elliptical hosts are gas-rich and therefore can- 
not be normal ellipticals. Our findings prove that not only 
elliptical hosts contain more gas, but that this gas has a 
significant level of ionization compared to quiescent ellipti- 
cals. However, the available data do not allow to determine 
unambiguously the ionization source. For instance, the ion- 
ization might be due either to the QSO radiation or to shock 
waves during a merger event. 

5.3 Magnitude- Asymmetry relations 

If mergers and/or gravitational interactions are to trigger 
or enhance QSO activity, we might expect a correlation 
between the QSO magnitude and the degree of asymme- 
try of the system. The relation is shown in Fig. |9] for el- 
lipticals and spirals separately because they have a differ- 
ent behaviour. Indeed, for the ellipticals, a reliable cor- 
relation {p — 0.54) exists between the absolute magni- 
tude and the asymmetry coefficient (solid line in Fig. |9]): 
MviQSO) = (-1.10 ± 0.36)a(V) ± (-23.22 ± 0.26). It is 
even more robust (p = 0.68, dashed line) if the three higher 
asymmetry systems are removed from the fit, with a relation 

MviQSO) = (-2.51 ± 0.61)a(V) ± (-22.67 ± 0.30). (5) 

A similar relation is not found for the host magnitude. Con- 
versely, the asymmetry of spiral hosts does not correlate with 
QSO magnitude, but does correlate with host magnitude as 
Mv{Host) = (-1.09 ± 0.34)aCl/) ± (-21.89 ± 0.24), with 
p = 0.71. 

For ellipticals, in the low asymmetry range [aiV) < 
0.7), asymmetry is mainly due to shifts between the posi- 
tion of the QSO (cqso) and the center of luminosity of the 
host (cl). To illustrate this statement, we plot on Fig. [10] 
the asymmetry according to cqso {o-qso{V), left plot) and 
the asymmetry according to cl (aHost(V^), right plot) versus 
the distance between both centers cqso and cl. While the 
distance is strongly correlated to aQso{V), it is not so much 
the case for anost{V), proving that, for elliptical hosts, the 
asymmetry is due to relatively symmetrical galaxies which 
are not centered on the nucleus position. In the framework 
of galaxy evolution via mergers, it is expected that ellipti- 
cals are created during major mergers involving galaxies of 
similar masses. Such mergers are also likely to trigger the 
activity of the nucleus. Thus, as the merger of the galaxies 
is faster than the merger of the black holes first orbiting 
around each other, it would not be surprizing to detect off- 
centered activity in apparently relaxed galaxies, where the 
two SBMH are still sufficiently far from each other. Alterna- 
tively, an off-centered activity might also be due to a black 
hole recoil with its accretion dis c after the fusion of a black 
holes binary. Simulations from IVolonteri fc Madaul (|2008l ) 
predict the existence and observability of such off-centered 
activity. 

For spirals, the observed correlation with a{V) concerns 
Mv{Host) instead of My (QSO). Tidal disturbances in spi- 
rals are believed to arise either from minor mergers between 
dif ferent mass g alaxies (4 : 1 to 10 : 1 or higher mass ratios, 
see lBournaud et al. (2004)) or in young mergers which have 
not evolved during enough time to reach their final state. 
The lack of correlation found between asymmetry and QSO 
magnitude supports the idea that young or minor merg- 
ers have not infiuenced the activity of the central nucleus 
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Figure 9. Asymmetry is plotted against QSO and host absolute 
magnitude in the V filter for ellipticals and spirals separately. 
Open symbols indicate the presence of interactions or mergers. 
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Figure 10. Left (resp. right): distance in kpc between the center 
of luminosity of the galaxy and the QSO position against the 
asymmetry coefficient computed on the QSO position (resp. the 
center of luminosity of the host), for elliptical galaxies only. 



yet. This idea match es well with the concluding remark in 
iReichard etHI (|2009h . who suggest that the period of black 
hole growth may be preferencially associated with the end 
stages of minor mergers. On the other hand, morphological 
disturbances are already visible and contribute to raise the 
asymmetry in proportion to their luminosity, giving rise to 
the observed correlation. 



6 CONCLUSIONS 

We have studied a sample of 69 QSOs, amongst which 60 
have been observed with the ESO NTT and SUSI2. Our 
QSO-host separation technique, based on the MCS decon- 
volution method, proves well adapted for that kind of data. 
A morphology classification method based on the asymme- 
try of the system has enabled us to compare the properties 
of mergers and non interacting hosts. More precisely, high 
asymmetry systems tend to have a higher degree of ioniza- 
tion and a bluer continuum than low asymmetry systems. 
Another interesting finding is that the QSO hosts tend to 
contain on average more ionized gas than quiescent galaxies. 
This trend is especially strong for ellipticals. It is consistent 
with the scenario in which AGN host galaxies lie in a transi- 
tion phase in betwee n the red and blue m odes in the colour- 
magnitude diagram (jMartin et al.ll2007h . as these emissions 
will make the elliptical QSO hosts appear bluer than their 
quiescent counterparts. Nevertheless, our data do not allow 
to firmly determine if the ionization source is the QSO itself 
or shocks induced during merger events. 
The relation between the QSO and host magnitudes is quite 
weak, and mainly found in ellipticals. The correlation ob- 
served is probably only a remnant of the more fundamental 
black hole-spheroid relation. We also find that high asym- 
metry systems cover a wider range of host magnitudes at a 
given QSO magnitude, which is interpreted either as lumi- 
nosity excess due to galactic fusion or to particularly pow- 
erful nuclear activity. 

At this stage, we cannot firmly determine unequivocally 
which process is responsible for QSO triggering. Our study 
enlightens links between interaction processes, QSO activity 
and gas ionization, but it would be risky to assess unequiv- 
ocally any hierarchical structure between those process. In- 
deed, our study does not allow to detect explicitely gas in- 
falling to the central regions due to the interaction process. 
However, we find that the nuclear activity during a merg- 
ing process is different for ellipticals and spirals, reinforcing 
the current belief that they correspond to different stages 
of evolution in a QSO lifetime, which fits a merger-driven 
evolution scheme. 

A few ellipticals have an off-centered activity which might 
be due to a black hole recoil or to a merger induced activity 
between similar mass galaxies. Moreover, some particular 
cases are reported for the first time and described in more 
detail in the Appendix. They seem to have experienced a 
wide variety of merging scenarios, and certainly deserve fur- 
ther investigations. Finally, we report the first image of the 
underlying host of PDS456, and discuss the possibility that 
this system may contain a double QSO. 

As future works, the study of a sample of inactive galax- 
ies in similar redshift and brightness ranges with similar spa- 
tial resolution and filters will help clarifying several open 
questions, such as the frequency of mergers or the spread 
in colours compared to templates. A detailed comparison of 
the newly introduced asymmetry coefficient to other existing 
asymmetry measurements is also planned. 
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Figure Al. Observations of QSO 1307 + 085 in both filters. The 
top row (resp. bottom) shows, from left to right, an observation, 
the deconvolved image (with the special features discussed encir- 
cled), and the residual map for the V#812 filter (resp. WB#665). 
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APPENDIX A: PECULIAR CASES 

One of the goals of the present study is to find candi- 
date hosts which would represent a particular stage of evo- 
lution in a QSO's lifetime (for example, s imilar to the 
now f a mous HE0450 - 2958 case stu di ed in Magain et al.l 
ll2005l): iLetawe et al.1 (|2008bl . l2009h : lElbaz et al.1 l)2009l l: 
Ijahnke et al.l ( 20091 )). in preparation for higher resolution 
observations which would go deeper into the central regions 
and would allow to better assess the nature of QSO-host in- 
teractions. In the present Appendix, we review some special 
cases, which are essentially outliers from the main trends 
described in the previous section. Note that all images have 
a zero position angle (North is up, and East is to the right), 
and all indicated distances are projected distances. 



Al 1307+085 

First of all, the host of this QSO lies amongst the very few 
systems which have substantially different asymmetry co- 
efficients in the two filters {a{V) = 1.29, a{WB) = 0.28). 
Secondly, the very blue QSO colour Am{QSO) = —0.5 sug- 
gests that the nucleus has strong emission lines compared 
to the continuum. Moreover, the comparison of the decon- 
volved images in both filters (Fig. lAip reveals two zones 
(encircled in the figure) in the host that show a prominent 
emission in the V-filter not seen in the continuum, which 
we thus associate to strong emission lines. Those elements 
put the host in the subclass of asymmetric ellipticals with 
substantial gas ionization, where the strong emission zones 
increase the asymmetry coefficient and the magnitude in the 
V-filter, whereas the stellar component looks like a rather 
regular elliptical. An emission zone appears in both filters 
at ~ 22.5kpc West from the nucleus, close to the left edge 
of Fig. lAll However we see no obvious trace of interaction 
between this feature and the host. 



Observation 



Q0022-2044 
V#812 z-0.170 



Deconvolution 
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Figure A2. Observations of QSO Q0022 - 2044 in both filters. 
The top row (resp. bottom) shows, from left to right, an observa- 
tion, the deconvolved image, and the residual map for the V#812 
filter (resp. WB#665) 



A2 Q0022-2044 

Q0022-2040 has also substantially different asymmetries in 
both filters {a{V) = 0.73, a{WB) = 1.45). The high asym- 
metry in the WB#665 filter is due to a shift of 2 pixels = 0.9 
kpc between the center of luminosity of the galaxy and the 
position of the nucleus, as shown in Fig. lA2l This shift, along 
with the very blue host colour Am{Host) — —1.94, indicate 
that Q0022-2044 is a typical example of the class of ellipti- 
cals with strong gas emission lines harbouring off-centered 
activity already discussed in Section 5.3. Let us also men- 
tion that the QSO is particularly bright in the WB#665 
filter (Am{QSO) — 0.41), suggesting a very red continuum 
and weak emission lines. 



A3 PDS456 

PDS456 is known as the most powerful QSO in the local 
Universe {z < 0.3). It is t hought to be a radio-quiet ana- 
logue of the famous 3C273 l|Schmidtlll963h . Spectra have al- 
ready been obtained in basically all wavelength ranges and 
K-band, radio and CO (l-O) images are also available (see 
lO'Brien et al] (|2005l) or I Yun et all (2004) for a complete re- 
view and references). It was suggested that it is in a crit- 
ical transition phase between an Ultra Luminous Infra Red 
Galaxy (ULIRG) and a QSO (|Sanders fc Mirabel [l99^ . as 
it shows an optical spectrum dominated by broad emission 
lines, large IR and X-ray luminosities, as well as a large 
dust /cold gas content. Moreover, UV and X-ray spectra re- 
veal the presence of decelerating cooling outflows probably 
driven by radiation or magnetic field. The deconvolved im- 
age (Fig. IA3|I shows two compact sources indicated by El 
and E2, situated a sim2 ercseconds S-W of the QSO, and 
compatib le with th e three blended sources of the K-band 
image m I Yun et al.l \2004 ). Deconvolution allows to disen- 
tangle the underlying host, which seems elliptically shaped. 
It is, to our knowkedge, the first image of the galaxy host- 
ing this QSO in the visible. Magnitudes and asymmetries are 
given in Table [3] for the whole system (El-|-E2-|-Host), and 
decomposition in the different components is given in Table 
\M\ The difference in magnitude between the host and the 
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m{QSO) 


m{Host) 


m{El) 


m{E2) 


V 


13.42 


17.56 


19.35 


19.27 


WB 


13.10 


16.43 


18.53 


18.34 




m(Sl) m{S2) 








V 


14.17 14.17 


17.56 


19.35 


19.27 


WB 


13.31 14.99 


17.10 


18.73 


18.54 



Table Al. Magnitudes of the different components of the decon- 
volution of PDS456 with one (top) or two (bottom) point sources 
at the center of the host. 



QSO (—4.14) is the highest of the sample, which indicates a 
particularly strong activity of the nucleus. 
Let us now focus on the residuals map shown in Fig. IA3l for 
the V filter. It has a structure typical of the presence of two 
blended point sources with similar intensities. Indeed, if a 
point source is to be fit where in fact there are two closely 
blended ones, the two peaks will not be correctly fit and this 
will result in a double hole in the residual map (model minus 
observation). If such a feature is real, and does not corre- 
spond to a PSF mismatch on that particular observation, 
we expect to find it in each exposure. Figure IA4I shows an 
average of the four residuals resulting from the simultaneous 
deconvolution of 4 observations. The double-hole structure 
is clearly present. We tested the hypothesis of a double nu- 
cleus by adding in the deconvolution process a second point 
source, along with the usual diffuse background. The result 
is displayed on Fig. lASI In the V filter, the two point sources, 
of comparable magnitude My ~ —27.45 are only separated 
by O.Slkpc and the host nearly disappears, with only a faint 
tail starting from the nucleus. For the WB filter, the separa- 
tion is only 0.22kpc, and a faint host, whose morphology can 
hardly be determined because of its faintness and its small 
angular size, is detected. Magnitudes in both filters for the 
two point sources fit are given in the bottom part of Table 

Let us now examine arguments in favour and against 
the double active nucleus hypothesis. The very good quality 
of the deconvolution with two point sources makes clear that 
there is at least a very compact source in both filters next 
to the "main" nucleus. Moreover, for the one point source 
fit, the value Am{Host) = 1.13 is the highest of our sample 
and surpasses the typical elliptical value, indicating an ab- 
normally red continuum. This strong compact continuum is 
better fitted with a point source, in coherence with an hypo- 
thetical AGN activity. The presence of two separate nuclei 
might explain its high luminosity, and is not inconsistent 
with the accepted scenario of a transitory system between 
ULIRG and QSO via a merging process. It could also ac- 
coun t for the except ionally broad abs orption (jO'Brien et al.l 
I2005D and emission (|Yun et al.ll2003 ) features observed. 
Conversely, if both sources were to correspond to AGN activ- 
ity, we would expect a common center in both filters, which 
is not really the case {Sc ~ 0.32kpc). Also, their magnitudes 
behave differently in each filter (see Table IXT1) . their relative 
intensities being much more different in the WB filter than 
in the V filter. Moreover, the value Am{Host) — 1.13 might 
be interpreted as due to a strong reddening of the host, 
which would match the ULIRG hypothesis. Thus, from our 
observations, it is risky to assess their true nature unequiv- 
ocally. 




Figure A3. Observations of QSO PDS456 in both filters. The 
top row (resp. bottom) shows, from left to right, an observation, 
the deconvolved image, and the residual map for the V#812 filter 
(resp. WB#665). 




Figure A4. Average residuals from the simultaneous deconvolu- 
tion of 4 observations of the QSO PDS456 in the V filter. The 
two holes in the residuals suggest the presence of a second point 
source. 

All in all, it is clear that PDS456 is an exceptional object, 
but maybe it is exceptional in a presently unexpected way. 
Deep high resolution optical imaging, or 3D integral field 
spectro scopy, processed with a similar method as the one 
used in lLetawe et al.l (|2008bl ). might help clarifying its sta- 
tus. 

A4 PG 1012+008 

From Fig. EH it is clear that PG1012 -I- 008 is involved in a 
merger containing at least 2 galaxies separated by ~ lOkpc, 
and a third one 20kpc N-W from the QSO, which might 
also be interacting gravitationally. The best fit is obtained 
by using 3 point sources in the deco nvolution, located at 
each galactic center position. However. iBahcall et al.l (Il997t) 
have already observed this system with HST and WFPC2, 
and no other point source was found. This difference may 
be explained in the following way. Given the difference in 
resolution between both observations (0.1" /pxl for WFPC2 
compared to 0.161"/pxl for SUSI2), a very compact and in- 
tense stellar emission just resolved by WFPC2 might look 
unresolved in our SUSI2 observations. At the QSO redshift, 
the range of sizes in which an emission zone would be re- 
solved with WFPC2 but not with SUSI2 is 0.2 - O.Skpc. 
We analysed the morphology of those two putative point- 
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Figure A5. Deconvolved image in both filters of PDS456, with 
two point sources very near from each other. The average resid- 
uals are considerably improved, without any significant structure 
remaining. 
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Figure A6. Deconvolution process for PG1012 + 008, where 3 
point sources are included, which don't necessarily correspond to 
any AGN-activity. 



like sources using the WFPC2 archives. Figure IA7I reveals 
that the point-like source nearest from the QSO looks like a 
compact emission zone whose FWHM is 4 about pixels. The 
QSO FWHM in that image being only 2 pixels, this emission 
is clearly not point-like. However, this 4 pxl width converts 
in SUSI2 pixels to 2.5 pxl FWHM, which is very close to 
the resolution fixed in the deconvolution process (2 pixels 
FWHM), making it look very much like a point-source at 
this resolution. Moreover, the flux ratio between the QSO 
and the compact emission is the same (~ 0.045) as the flux 
ratio between the QSO and the point source in the SUSI2 
images, reinforcing the hypothesis that the compact emis- 
sion in WFPC2 and the second point source in SUSI2 are the 
same object. The same arguments hold for the third point 
source. Another indication that the additional compact re- 
gions are not related to any AGN activity is that the colours 
computed for the nearest galaxy center and the other one 
are Am = 0.57 and 0.62. Those values are higher than for 
any QSO in our the whole sample, revealing an important 
contribution of the continuum compared to emission lines. 
This fact tends to favour the view that the 2 added point 
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Figure A7. PG1012 -|- 008 observed with HST/WFPC2 and a 
zoom on the companion galaxy containing a compact, but not 
point-like, emission zone. 




WB#E65 



Figure A8. Deconvolution of 115H-117, with 3 point sources. 
The S-W one is probably a foreground star, whereas the two oth- 
ers are hosted by interacting galaxies. 



sources do not correspond to AGN activity, but rather to 
compact luminous stellar regions, characterized by a promi- 
nent continuum. 



A5 1151 + 117 

Apart from the central QSO, the observation of 1151-1-117 
(Fig. IA8|I reveals the presence of two emitting regions 
20.6kpc S-W, and 11.5kpc N-E of the QSO. Both need 
a point source component to achieve good residuals. De- 
convolution reveals that the S-E one is well flt by a sin- 
gle point source, suggesting it might be a foreground star. 
On the other hand, the N-E one needs a smooth back- 
ground to achieve a good flt. This point source is quite weak 
(My = —19.32), and has a dominant continuum emission 
Am = 0.33. Similarly to PG1012+008, this high value of 
Am might indicate that the point-like emission is due to a 
compact stellar region, even if it is still within the range of 
reasonable value for an active nucleus. 

The background needed around the second point source 
shows clear signs of gravitational interactions with the host, 
especially on the WB#665 deconvolved image, where the 
galaxies seem to share a tidal tail, suggesting a conflguration 
similar to the well known M51 (jSalo fc LaurikainenI l2000l : 
iToomre fc Too mrc 1972), where two galaxies of mass ratio 
3:1, viewed nearly face-on, probably encountered for the 
flrst time 300 ± 100 Myr ago. Two additional faint emissions 
are also present in both fllters (N-W and N), and may also 
be linked to the merging process. 
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Figure A9. Deconvolution process for HE1202 — 0501. The very 
atypical morphology may be created by the collision of two galax- 
ies of different disc mass. The rather poor quality of the residuals 
is likely due to significant short-scale structures near the center. 



A6 HE1202-0501 

The host of HE1202-0501 (Fig. \M} is very atypical, with a 
45 kpc nearly straight tail extending from the center towards 
the S, and containing a few minor substructures. Another 
much weaker tail of ~ 27 kpc length extends from the center 
to the N-W. Such a special configuration is very similar to 
the shape of a 1.2 — 1.3 Gyr old major m erger of nearly equal 
mass disc galaxies (|SpringeI et al.l (|2005l ). Fig. 7). These sim- 
ulations also show very complex intense features next to the 
center. If such structures are prensent in HE1202-0501, they 
might not be resolved by SUSI2, but could substantially 
contribute to the observed total central nucleus luminosity, 
leading to PSF mismatches and thus decreasing the quality 
of the deconvolution around the center, as can be seen in 
the residuals shown in Fig. \M\ HE1202 - 0501 also looks 
very similar to the fa mous "Mice" merger, NGC4746, which 
is also well rnodelle d (|Toomre fc Toomrelll972l: iMihos et al.l 
1 19931 : iBarne^ l2004l ) by a two identical mass disc galaxies 
merger event, where one of the two discs is viewed edge-on 
and mimics an elongated bar. They are a few major differ- 
ences between NGC4746 and HE1202-0501. In "The Mice", 
the galactic centers are separated by 20kpc, and both tails 
have approximately the same length and luminosity, whereas 
HE1202-0501 shows only one nucleus, and different size tails. 
Different tail size and luminosities are suggestive of differ- 
ent mass galaxies. We thus conclude that IIE1202-0501 is 
most probably a merger between two disc galaxies of differ- 
ent masses containing a sufficient amount of gas to induce 
AGN activity, one of the two discs being viewed edge-on. 
The poor residuals near the center may result from intense 
star formation activity around the nucleus. 

A7 1023-014 

The QSO 1023 — 014 is involved in a major merger involv- 
ing two galaxies of comparable size. The best deconvolu- 
tion is achieved with 3 point sources, labelled A, B and 
C on Fig. lAlOl While A and B have are very similar in 
both filters {mv{A) = 19.54, mv{B) = 19.67, Am{A) = 
—0.06, Am{B) — —0.01), C has a quite different behaviour, 
with mv{C) = 20.55 and Am(C) — 0.51, which deviates 
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Figure AlO. Deconvolution process for the major merger of the 
QSO 1023 — 014. 3 point sources are necessary to achieve satis- 
factory residuals. 

strongly from the mean Am{QSO). As in the case of PG 
1012 -I- 008, it indicates that the source C might not be re- 
lated to AGN activity but rather to strong unresolved stellar 
emission. Moreover, this system lies significantly below the 
Mv {QSO)- My (Host) relation, having an underluminous 
AGN compared to its host. Mv{QSO) - Mv{Host) = 1.27 
is the highest value of the whole sample, the nuclear magni- 
tude being indeed below the standard delimitation between 
QSO and Seyfert. All in all, and given the fact that the 
QSO and total host are more powerful in the WB#655 fil- 
ter (the one containing no emission line), 1023 — 014 could 
be a rather dry merger, with no sufficient amount of gas to 
create powerful AGN activity, but possibly one or two low 
luminosity AGNs. Further studies, including spectroscopy of 
the nucleus, would be necessary to test this hypothesis. 

A8 HE1211-1905 

The QSO HE1211-1905 is surrounded by a spiral galaxy 
which shows clear signs of gravitational interaction via ex- 
tended tails N and E of the nucleus, corresponding to the 
Tl and T2 boxes in Fig. lAllI Given the significant area of 
those tails (~ lOkpc^), it makes sense to check differences 
in magnitudes between both filters in those regions in order 
to specify their nature. 

The analysis reveals that Tl and T2 are very similar, as 
they both have Am — —0.05. This negative value is much 
lower than the value computed for the whole host (0.5). It 
indicates that ionized gas is more prominent in the faraway 
regions Tl and T2 than in the host's central regions. The tail 
Tl goes even beyond the extracted image. Fig. IXT2I shows 
a larger area around the QSO, where the size of Tl can be 
estimated to ~ 40 kpc. The cause for such disturbances is 
unclear because HE1211-1905 appears quite lonely, appart 
from a tiny galaxy, probably elliptical, lying 35 kpc from it. 
Let us mention that there are also a few spots in the over- 
all direction of that companion galaxy, which may also be 
linked to the system. Another possibility is that the system 
is a late-stage merger, where the two galaxies nuclei have 
rapidly merged, leaving a host not yet relaxed. 
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Figure All. Deconvolution process for the strongly disturbed 
host of the QSO HE1211-1905. 




Figure A12. Surrounding of HE1211 - 1905. The tidal tail of 
~ 40kpc is surprisingly long. A companion galaxy and a few spots 
which arc present on the same side may be linked to the host 
disturbances. 



